Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 7 May 2009 (MN KTbX style file v2.2) 



Circumplanetary disc properties obtained from radiation 
hydrodynamical simulations of gas accretion by protoplanets 

Ben A. Ayliffe* and Matthew R. Batet 

School of Physics, University of Exeter, Stocker Road, Exeter EX4 4QL 



7 May 2009 



ABSTRACT 

We investigate the properties of circumplanetary discs formed in three-dimensional, self- 
gravitating radiation hydrodynamical models of gas accretion by protoplanets. We determine 
disc sizes, scaleheights, and density and temperature profiles for different protoplanet masses, 
in solar nebulae of differing grain opacities. 

We find that the analytical prediction of circumplanetary disc radii in an evacuated gap 
(^Hill/3) from Quillen & Trilling | ( 1998 ) yields a good estimate for discs formed by high mass 
protoplanets. The radial density profiles of the circumplanetary discs may be described by 
power-laws between r"^ and r We find no evidence for the ring -like density enhancements 
that have been found in some previous models of circumplanetary discs. Temperature profiles 
follow a ~ power-law regardless of protoplanet mass or nebula grain opacity. The discs 
invariably have large scaleheights {Hjr > 0.2), making them thick in comparison with their 
encompassing circumstellar discs, and they show no flaring. 

Key words: planets and satellites: formation - accretion, accretion discs - hydrodynamics - 
radiative transfer - (stars:) planetary systems: formation - methods: numerical 



1 INTRODUCTION 

The giant planets of our solar system all have multiple satellites, di- 
vided into two classes; the regular and irregular satellites. The regu- 
lar satellites have low eccentricity orbits with low inclinations rela- 
tive to the planet's equatorial plane ( jGalileij 1610||1989{ |Mosqueira| 
|& Estrada|2003l l, and possess prograde orbits. These characteris- 
tics are suggestive of satellite formation in a disc of material ro- 
tating about the planet. Contrastingly, the irregular satellites have 
eccentric, highly inclined orbits, and show no strong bias between 
prograde or retrograde orbits ( |Grav et al.|2003| l. These character- 
istics suggest that the irregular satellites are captured by their host 
planets, rather than formed in situ. However, this capture process 
requires a mechanism by which to lose energy, for example gas 
drag as suggested by |Pollack, Bums & Tauberj ( |1979[ >. This was 



Cuk & Bums 



1 2004b, who modify the 



more recently revisited by 
original scenario to include gas drag from a circumplanetary disc 
as opposed to the spherical envelope assumed by Polla ck et al.| 



Cuk & Burns find that drag due to the circumplanetary disc gas is 
sufficient to enable the capture of bodies up to 100km in radius, 
which would otherwise form only temporary associations. There- 
fore circumplanetary discs may be instrumental in the formation of 
satellite systems comprising both regular and irregular satellites. 
To form regular satellites a disc must exist around a proto- 
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planet in which agglomeration processes can occur. [Canup & Ward| 
( ,2002 l discuss the uncertainty in the nature of such a disc, whether 
it is an accretion circumplanetary disc or a spin-out disc. The for- 
mer is created when a protoplanet begins to contract whilst still 
embedded in the circumstellar disc. When the protoplanet attains 
a sufficiently small radius it becomes impossible for the high an- 
gular momentum gas from the circumstellar disc to fall directly 
onto its surface. Instead it forms an accretion disc through which 
angular momentum can be dispersed before the gas infalls to the 
protoplanet's atmosphere. A spin-out disc might form when a pro- 
toplanet contracts after the circumstellar disc has dispersed. In this 
scenario, as the protoplanet contracts, (as there is no ambient envi- 
ronment into which the angular momentum may be disseminated) 
so its rotation rate increases. Upon becoming rotationally unstable, 
the most rapidly rotating material which occurs at equatorial lati- 
tudes, is flung out to form a spin-out disc. 

One-dimensional models suggest that a Jupiter mass proto- 
planet will contract to less than 1/lOOth of its Hill radius in un- 
der IMyr (Papaloizou & Nelson 20051. This would be within the 
circumstellar disc lifetime, allowing the protoplanet to continue to 
accrete material via an accretion disc. This suggests that it is an ac- 
cretion disc in which satellites form rather than in a spin-out disc. 
Indeed, previous two-dimensional (Lubo w~Seibert & Artymow-| 
icz 1999; D'Angelo, Henning & Kley 2002) and three-dimensional 
(|Bate et al.H2003( |D' Angelo, Kley & Hennmgl[2003| [Fouchet &| 
|Mayer||2008[ |Ayliffe & Bate||2009[ > hydrodynamical models of 
planet formation have shown the existence of circumplanetary ac- 
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cretion discs during giant planet formation. The two-dimensional 
models that develop circumplanetary discs have shown strong 
shocks within these discs l [Lubow et al.|1999||D' Angelo et al.|2002^ 
that rapidly deplete them of material, which accretes on to the pro- 
toplanet. |Bate et aL]j2003^ showed that in three-dimensional mod- 
els these shocks appear to be much weaker than in two dimensions, 
suggesting a slower rate of circumplanetary disc depletion. Such 
weakened shocks also appear in recent three-dimensional calcula- 
tions by |Machida| ( (2008| . Following dispersal of the circumstellar 
disc, the circumplanetary disc is no longer replenished with gas and 
solids. The residual circumplanetary disc must presumably survive 
for a considerable period to allow satellites to form from the volatile 
compounds that will only solidify at low temperatures which can- 
not be achieved in an actively accreting disc. A long disc lifetime 
also allows time for the constituent materials (ice, rock, etc.) to ag- 
glomerate, and to do so slowly enough that the collisional heating 
does not sublimate the ice ( jCanup & Ward|2008^ . 

In this paper, we present and discuss the properties of cir- 
cumplanetary discs formed in three-dimensional radiation hydro- 
dynamical calculations of giant planet growth. These models su- 
persede the previous three-dimensional models mentioned above in 
that they include radiative transfer Therefore, they provide a more 
accurate description of the early phase of circumplanetary discs, in 
which satellites are believed to form. We examine the radial extents, 
density profiles, temperatures, and scaleheights of discs around dif- 
ferent mass protoplanets, in nebulae with different grain opacities. 
In Section 2, we describe our computational method. Section 3 
presents the results from our self-gravitating radiation hydrody- 
namical simulations, while in Section 4 we discuss our results in 
the context of previous work. Our conclusions are given in Section 
5. 



2 COMPUTATIONAL METHOD 

The calculations described herein have been performed using a 
three-dimensional smoothed particle hydrodynamics (SPH) code. 
This SPH code has its origins in a version first developed by |Benz| 
( [19901 [BenzeTaLlfiggO) but it has under gone substantial modifi- 
cation in the subsequent years. Energy and entropy are conserved 
to timestepping accuracy by use of the variable smoothing length 
formalism of |Springel & Hemquist| p002| > and |Monaghan j2002^ 
with our specific implementation being described in Price & Bate| 
( |2007^ . Gravitational forces are calculated and neighbouring par- 
ticles are found using a binary tree. Radiative transfer is modelled 
in the flux-limited diffusion approximation using the method devel- 
oped by|Whitehouse, Bate & Monaghan|(2005[ > and |Whitehouse &] 
|Bate| j2006 l|. Integration of the SPH equations is achieved using a 
second-order Runge-Kutta-Fehlberg integrator with particles hav- 
ing individual timesteps (Bate, Bonnell & Price 1995). The stan- 
dard implementation of artificial viscosity is used (Monaghan &^ 
|Gingold| 1 9"83| |Monaghan| 1 992^ with the parameters a,, = I and j8„ 
= 2. The code has been parallelised by M. Bate using OpenMP. 

The calculations used in this paper have been performed in the 
same way, and in some instances are the same as calculations per- 
formed in |Ayliffe & Bate| ( [2009) (henceforth Paper I). It is therefore 
only in brief that the details of the method are outlined below, with 
more details available in Paper I. 



2.1 Model setup 

Our calculations are performed in the reference frame of the planet, 
orbiting a star of mass M, at radius with an angular speed given 

by Q.p = ^GM,/r^, neglecting the mass of the planet. The stellar 
mass is taken as 1 Mq and the orbital radius of the planet to be 5.2 
AU. We model only a small section of the protoplanetary disc cen- 
tred on the protoplanet. Our standard section size is r = 1 ± 0. 15 fp 
(5.2 ± 0.78 AU), and (p = ±0.15 radians. The initial radial tem- 
perature profile for the disc is oc r^' , which leads to a constant 
ratio of disc scaleheight with radius of H/r = 0.05. The initial sur- 
face density of the disc has as E oc r"'/^ profile, with a value of 
75 g cm"^ at the planet's orbital radius. 



2.2 Boundary conditions 

The distribution of particles within the disc section is initially that 
of an unperturbed disc, with Keplerian velocities. Inward and out- 
ward of the planet's orbital radius the particles flow in opposite 
directions in the reference frame of the planet, and therefore leave 
the modelled section by opposite azimuthal boundaries. Particles 
are injected to replicate the natural flow expected of a complete 
disc. The distribution of the injected particles is derived from the 
isothermal calculations of [Bate et al . (20031, interpolated to the 
pertinent protoplanet mass where necessary. This acts to replicate 
the expected gap in the vicinity of the protoplanet that would form 
were the entire disc modelled. The fact that the ZEUS calcula- 
tions are vertically isothermal leads to a small mismatch in the 
vertical temperature structure at the boundaries for the radiation 
hydrodynamical calculations. These differences are small, particu- 
larly when compared with temperature changes that occur near the 
planet due to gas accretion; see Paper I for more details. 

Ghost particles are employed along the boundaries of the disc 
section to provide the pressure forces that would come from the gas 
which is not modelled outside of the boundaries. 

A boundary region above and below the disc, located at a 
height equivalent to an a optical depth of 1 (t » I) from infin- 
ity, allows the radiative transfer scheme to transport energy from 
the bulk of the disc to the boundary as though radiating to infinity. 
SPH particles within the boundary regions are evolved normally, 
except that they interact with the SPH particles in the bulk of the 
disc without their temperatures being affected. 



2.3 Equation of state 

We use an ideal gas equation of state for our radiation hydrodynam- 
ical calculations; p = pTgRg/jj, where Rg is the gas constant, p is 
the density, Tg is the gas temperature, and /i is the mean molecular 
mass. The radiative transfer acts such that work and artificial vis- 
cosity (which includes both bulk and shear components) increase 
the thermal energy of the gas, and work done on the radiation field 
increases the radiative energy which can be transported via flux- 
limited diffusion. Energy transfer between the gas and radiation 
fields depends upon their relative temperatures, the gas density, and 
the gas opacity. For more details of the implementation see |White^ 
[house et al.| l |2005| l and |Whitehouse & Bate1j2006| l. 



2.4 Opacity treatment 



We use interpolation through the opacity tables of [Pollack et al.| 
( |1985 1| to provide the interstellar grain opacities (IGO) for so- 
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lar metallicity molecular gas, whilst the gas opacities, which be- 
come important at high temperatures, are obtained from the ta- 
bles of Alexander (1975 ) (the IVa King model); for further details 
see [Whitehouse & Bate (20061. The agglomeration of interstellar 
grains into larger grains, or the sublimation of small grains, may 
lead to grain opacities that are lower than the IGO ( Podolak'2003 1. 
We mimic these elfects by dividing the IGO values of Pol lack et al.| 
(jT985} by factors 10, 100, or 1000 (details in Paper I). 



2.5 Protoplanets 

Our protoplanets are modelled by a gravitational potential, and a 
surface potential that yields an opposing force upon gas within one 
protoplanet radius of the protoplanet surface. The combination of 
the gravitational and surface forces takes the form of a modification 
to the usual gravitational force as 



F, = -- 



GM„ 



1 - 



2R„ - r 



(1) 



for r < 2 Rp where r is the radius from the centre of the proto- 
planet, Rp is the radius of the protoplanet, and Mp is the mass of 
the protoplanet. This equation yields zero net force between a par- 
ticle and the protoplanet at the surface radius Rp, whilst inside of 
the protoplanet's radius the force is outwards and increases rapidly 
with decreasing radius. 

The majority of our calculations use protoplanets with radii 
equivalent to 1% of their respective Hill radii (at 5.2 AU), but in 
one instance we model a Jupiter mass protoplanet with a radius 
equivalent to present day Jupiter (7.15 x 10' cm). A smooth start to 
the radiation hydrodynamical calculations is given by exponentially 
shrinking the protoplanet radii from an initial value of 0.01 rp down 
to the desired radii during the first orbit of the protoplanet. The 
models settle to have steady temperature and density profiles after 
less than 2 orbits, with changes thereafter being very gradual (see 
Paper I, Figure 4). 



3 RESULTS 

We modelled protoplanets ranging in mass from 10-333 M0. 
Figure [T| illustrates the density distribution about 100, 166, and 
333 M0 protoplanets in both the x - y and x - z planes. We fo- 
cus here on the immediate vicinity of the protoplanet. The reader 
interested in the structure on larger scales such as the spiral arms 
and gap in the circumstellar disc should refer to Paper I. The dis- 
tributions about a 100 M0 protoplanet in a nebula with interstellar 
grain opacities (IGO) are shown in the top-left panels. The density 
distributions in the x - y and x - z slices are very similar; this case 
appears to be spherically symmetric, and the protoplanet is seen 
to possess no circumplanetary disc. However, the density distribu- 
tions about a 100 M0 protoplanet in a 1% IGO nebula, shown in 
the bottom-left panels of Figure[T| make evident the existence of a 
circumplanetary disc. 

To determine the size of the circumplanetary discs formed in 
our models we adopt a specific set of criteria. The radial edge of 
a circumplanetary disc is taken as the point of turnover in the spe- 
cific angular momentum of the disc, calculated from the material 
within 5° of the midplane. The top panels of Figure |2] show the 
specific angular momentum distributions of gas around different 
protoplanets, in nebulae of different grain opacities. Furthermore, 
the absolute value of the angular momentum at the turnover point 



can be no more than a factor 2.2 different from the expected Keple- 
rian specific angular momentum (denoted by straight dashed lines 
in Figure [2| at the given radius. The absolute difference in the spe- 
cific angular momentum is due to pressure support which reduces 
the velocity that the gas must have to maintain a particular orbital 
radius. 

Using the described criterion to determine whether a particular 
protoplanet has a disc, and the size of that disc, we found that only 
the 100, 166 and 333 protoplanets possess circumplanetary 
discs. Figure [3] illustrates the circumplanetary disc radii obtained 
for these protoplanets. Protoplanets of lower mass exert an insuffi- 
cient vertical gravity component to overcome the thermal pressure 
in their surrounding envelopes, thus they remain spherical. Just as 
was seen in Figure [T| the 100 M0 protoplanet in an interstellar 
grain opacity nebula is not found to have a circumplanetary disc. A 
simple analytic prediction of the approximate circumplanetary disc 
radii around accreting protoplanets was made by Quill en & Trilling | 
([1998); their derivation is as follows. Assuming that the protoplanet 
is embedded in a gap, as is the case for the large masses described 
here, the mass inflow into the Hill sphere will enter via the inner 
and outer Lagrange points. An estimate of the specific angular mo- 
mentum of the gas at these points can be made by assuming that 
the gas' transverse velocity component, v^, comes purely of the 
Keplerian shear. In the rotating frame of the star-planet system 



and the Lagrange points are approximately at 

r=rp± Rmu 

where the Hill radius is 



^Hill — i\ 



Thus, 



Mp 



1/3 



1 + ■ 



(2) 



(3) 



(4) 



(5) 



Applying an expansion to the term with the form (1 -I- x) ''^ 
equation[5]we obtain 



GM. 



(6) 



Thus, an expression for the specific angular momentum of gas at the 
Lagrange points relative to the protoplanet is j ^ RmV^'i. — ^mu^p- 
Assuming conservation of angular momentum, we estimate the 
centrifugal radius, R^, of the circumplanetary disc as 



]_ 



GMp 
Rc 



which can be solved to give 
Re ~ ^Hill/3. 



(7) 



(8) 



In all cases where circumplanetary discs are found, the disc radii 
do not differ greatly from this Rum/ 3 estimate for a protoplanet 
embedded in a gap (Figure|3j. It is interesting to note that reducing 
the nebula grain opacity may enable disc formation, but that beyond 
that, changing the opacity appears to have only a small impact on 
the resulting circumplanetary disc radius. 

In the second row of Figure |2] we show the midplane density 
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Figure 1. Midplane and vertical density profiles of the gas distributions surrounding protoplanets with masses of 100, 166, and 333 (left to right) and 
two different grain opacities. Top row: Density at the midplane in a nebula with interstellar grain opacity. Second row: Density in a vertical slice through 
the protoplanets, corresponding with row above. Third row: Density at the midplane in a nebula with 1% interstellar grain opacity. Bottom row: Density in a 
vertical slice through the protoplanets, con'esponding with row above. The spatial dimensions in units of the orbital radius (5.2 AU). 



profiles of each circumplanetary disc. Excluding the inner most re- 
gion where the density turns up due to gas piling upon the pro- 
toplanet, each disc density profile can be approximately described 
by a single power law. The exponent varies between -3/2 and -2 
for different mass protoplanets and for nebulae of different grain 
opacities. The change of exponent with protoplanet mass and neb- 



ula opacity does not show an illustrative trend in either regard. The 
surface density distributions of the circumplanetary discs are shown 
in the third row of Figure |2] and follow a similar downward trend 
with radius, but with a shallower falloff (~ r"''^ consistent with a 
constant H/r; see below). 

The bottom row in Figure |2] shows the temperature profiles 
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0.01 0.10 0.01 0.10 0.01 0.10 

Radius (Hill Radii) Radius (Hill Radii) Radius (Hill Radii) 

Figure 2. Specific angular momentum (top row), midplane density (second row), surface density (third row), and midplane temperature (bottom row) profiles 
of the gas distributions surrounding protoplanets with masses of 100, 166, and 333 Mgj (left to right) in nebulae of dift'erent opacities; IGO (solid line), 
10% IGO (dotted), 1% IGO (dashed), and 0.1% (dot-dashed). The straight dashed line in each panel in the top row indicates the expected specific angular 
momentum for material in Keplerian orbit about the protoplanet. The vertical dot-dash lines mark i?Hiii/3. 



© 0000 RAS, MNRAS 000, 000-000 



6 B.A. Ayliffe £r M.R. Bate 



0.6 



0.5 



0.4 



0.3 



0.2 



0.1 



0.0 



—I I I I ' ' ' ' ' 



_l I l_ 



10 100 
Opacity Divisor 



1000 



Figure 3. The radial extent of circumplanetary discs around 100 (dash- 
dots-dash line), 166 (dashed line), and 333 (solid line) protoplanets 
against varying grain opacity for the nebula material. 



calculated by azimuthally averaging the temperature of gas at the 
midplane of the circumplanetary discs. The temperatures are rel- 
atively high, ranging from 170K to 1600K. The lowest tempera- 
tures are achieved at reduced grain opacities. The fall off in tem- 
perature with radius for all opacities and protoplanet masses is ap- 
proximately described by a power-law, r"'''". For our Jupiter ra- 
dius protoplanet calculation the disc's inner edge is nearer to the 
protoplanet's centre, and thus within a steeper potential leading to 
a correspondingly higher temperature, up to ~ 4500K. The outer 
edge has a temperature of 240K, only slightly larger than obtained 
for the same calculation performed with a 1 % i?Hiii surface radius 
(Figure|4|. Thus the chosen size of the protoplanet has only a small 
elTect on conditions further out in the circumplanetary disc, and 
does not change our main conclusions. 

All of the circumplanetary discs that are produced in our mod- 
els are thick. The scaleheights (H/r) are shown in Figure[5]and can 
be seen to range between ~ 0.3 - 0.6. The scaleheights are con- 
stant with radius, showing no kind of flaring. Scaleheight tends to 
decrease with decreasing grain opacity due to the greater radiative 
cooling possible at lower opacities. This leads to cooler discs (as 
can be seen in Figure[2]( and therefore more flattened circumplane- 
tary discs (as can be seen in Figure[TJ. Only the lowest grain opacity 
does not follow this trend as the gas and radiation can decouple at 
such low opacities, reducing the effectiveness of radiative cooling. 
The scaleheights also tend to decrease with increasing protoplanet 
mass due to the larger vertical gravity component; once again, this 
effect is visible in the x - z slices in Figure [T| The scaleheight for 
the disc about our Jupiter radius protoplanet is broadly in line with 
the other 333 M0 protoplanet scaleheights, however it is slightly 
thinner than its 1% ifniu surface radius counterpart; H/r ^ 0.3 and 
H/r ^ 0.4 respectively. 

The X - z plane density plots in Figure[T|also include velocity 
vectors, to indicate the relative velocity of gas flow in the vicinity 
of the protoplanet. The gas velocities down on to the disc and pro- 
toplanet from high latitudes are considerably greater than the ve- 
locities through the disc. This velocity structure, with fast infall of 




0.0- 

Radius (Hi 



Figure 4. The same quantities as plotted in Figure |2] but for our 333 
protoplanet with a surface at a Jupiter radius instead of 1% Rn\u - The dot- 
dash line gives the results for the 1% /?Hiii case at an equivalent time for 
comparison; this is at an earlier time than the 1% /?Hiii case in Figure|2]as 
the small core calculations are slower, and thus less evolved. The nebula has 
1% interstellar grain opacities. 
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Figure 5. Scaleheights of circumplanetary discs around 100 (dash- 
dots-dash line), 166 (dashed line), and 333 (solid line) protoplanets 
against varying grain opacity for the nebula material. 




Angle from z-axis (radians) 

Figure 6. We investigate the latitudinal dependence of the mass inflow rate 
through the Hill sphere on to a 333 protoplanet embedded in a nebula 
with interstellar grain opacity. We plot the normalised cumulative distribu- 
tion of mass flux against latitude (the disc midplane is at nIT). The vast 
majority of the mass flows into the Hill sphere near the equator. 

gas from high latitudes, reflects what was found by |Klahr & Kley| 
l |200 6'l. However, the gas above and below the midplane is very 
diffuse, thus the mass flux in this direction contributes a very small 
fraction to the growing protoplanet. This can most clearly be seen 
in Figure [6] which illustrates that a majority of the mass enters the 
Hill sphere of the protoplanet at around 90° from the z-axis (i.e. 
near the x-y plane), feeding in to the circumplanetary disc. 



4 DISCUSSION 

The properties of circumplanetary discs around protoplanets are 
important in understanding the formation and capture of the satel- 
lite systems we see around giant planets in our solar system. 

4.1 Comparison with previous work 

A massive protoplanet that establishes a well defined gap in its par- 
ent circumstellar disc is expected to accrete material that flows into 
that gap via its Lagrange points. This gas in turn is expected to have 
low velocity relative to the protoplanet upon entering the gap due 
to its similar orbital radius with respect to the star. Using these as- 
sumptions |Qunien^|TOnmg||T998^ estimated a centrifugal radius 
of ~ i?Hiii/3; at this radius gas is expected to enter into orbit of 
the protoplanet. Gas must then evolve towards the protoplanet by 
angular momentum transport, forming a circumplanetary disc of 
material as it slowly infalls. The radial extents of the circumplane- 
tary discs in our models are in good agreement with this estimate, 
and markedly larger than those assumed in the models of |Canup"&| 
Ward (2002 ). Quillen & Trilling also give an estimate for the ex- 
pected disc radii about protoplanets that have not evacuated a disc 
gap. In these instances the gas flowing into the vicinity of a pro- 
toplanet has lower specific angular momentum, and the resulting 
circumplanetary disc is considerably smaller. However, we find no 
clear evidence of circumplanetary discs around protoplanets of less 
than 100 M0, instead we find near spherical envelopes. 

[D'Angelo, Herming & Kley| ( |2003^ performed two- 
dimensional hydrodynamical calculations to investigate cir- 
cumplanetary discs around protoplanets, including viscous heating 
and simplified radiative cooling. For a Jupiter mass protoplanet, 
they also find circumplanetary discs with outer edge radii (as 
defined by a turnover in specific angular momentum) of around 
/?Hiii/3, inline with |Quillen & Trilling| ( [T998 ) and our own results. 
The initial properties of our circumstellar disc (in which our 
planets are embedded) are similar to those of [D'Angelo et al.| 
(|2003j, with an identical surface density falloff, and similar initial 
temperatures and densities at the planet's orbital radius. The 
circumplanetary disc properties that develop in our calculations 
are also similar. For our Jupiter mass protoplanet, the circumstellar 
disc at standard opacity reaches a maximum temperature of 1600K 
at a 0.02/?Hiib whilst [D'Angelo et al.| find a maximum temperature 
of 1500K at the same radius. The surface density profiles in both 
sets of models are also similar. 

[Machida^ (^2008, 1 found that the circumplanetary discs formed 
in his calculations possessed a peak in surface density around 30 
protoplanet radii from the protoplanet's centre. Figure |4] shows the 
surface density and azimuthally averaged midplane density against 
increasing radius for our calculation of a Jupiter mass protoplanet 
with a Jupiter radius. The midplane density, and likewise the sur- 
face density, always decrease with increasing radius from the pro- 
toplanet out to the Hill radius with no sign of a ring-like density 
enhancement. This difference between Machida's results and our 
own is most likely a consequence of our differing treatments of the 
gas nearest the protoplanet. Machida uses a sink cell to remove gas 
at a radius approximately equal to a Jupiter radius, creating an evac- 
uated central region that leads to an inward pressure gradient near 
the protoplanet. Our more realistic approach is to model a surface 
upon which the gas piles up, increasing in density and temperature, 
and exerting an outward force on gas external to it. 

Our three-dimensional models show no evidence of strong 
shocks within the circumplanetary discs, in agreement with both 
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[Bate et al.H2003^ and |Machida|j2008^ , and in contrast to previous 
two-dimensional models. This again suggests that the strong shocks 
seen in two-dimensional models are produced by a lack of vertical 
spread in the colliding gas streams. A lack of strong shocks may 
be more conducive to the agglomeration of solids as the disc may 
be more quiescent. Shock driven migration of material through the 
discs may also be reduced, leading to longer lived discs. However, 
in contradiction to these previous three-dimensional calculations, 
we find very much thicker circumplanetary discs. This is caused 
by the trapping of heat in the optically thick circumplanetary disc; 
unable to radiate away the heat, the disc puffs up. A hotter disc 
may lead to faster accretion, opposing any enhanced lifetime gained 
from the weakened shocks. 

4.2 Implications for satellite formation 

The temperatures found in our circumplanetary discs are quite high, 
particularly when considering the large fraction of ice found in 
some of the Galilean moons. The inner most Galilean satellite, lo, 
is mostly composed of silicates and metals, Europa, mostly sili- 
cates with ~ 10% ice, whilst Ganymede and the outermost Galilean 
satellite, Callisto, both include some 50% ice ( |Schubertetal.|1986| i. 
The usual interpretation of the increasing ice fraction with orbital 
radius is of a radially decreasing temperature within the progenitor 
circumplanetary disc l |Lunine & Steve nson 1982). This temperature 
profile must reach a low enough temperature around Ganymede's 
orbit for large amounts of water ice to form and be incorporated 
in to the forming satellites ( jCanup & Ward|2002| ). However, only 
at the outermost reaches of our circumplanetary discs at reduced 
opacities does the temperature approach the ~ 160K expected for 
water condensation (at 5AU around a solar type star, |Stevenson| 
|1989[ l. This is at a radius some 17 times the present orbital radius 
of Ganymede. This may point towards orbital migration, though 
we have shown that there is a dearth of material at very large 
radii within a circumplanetary disc, which makes satellite forma- 
tion more challenging further from the protoplanet. It is also pos- 
sible that satellite formation occurs after significant gas accretion 
from the circumstellar disc has ceased, and the residual disc has 
cooled. Alternatively, perhaps ice must be delivered from the sur- 
rounding solar nebula in sufficiently large bodies such that it does 
not sublimate prior to satellite growth. 

It could also be that the moons were captured. However, the 
capture of the Galilean satellites is not supported by their low incli- 
nations and eccentricities, which are suggestive of formation within 
a circumplanetary disc. However some of Jupiter's satellites which 
have both low inclinations and eccentricities do appear to have 
formed in locations other than their present orbits, for example 
Amalthea ^Cooper et al. 2006). Amalthea orbits inside of lo, and 
yet has a density lower than water ( [Anderson et al.|[2005[ ), sug- 
gesting it formed elsewhere. However its low inclination and ec- 
centricity can be explained through its interactions with lo jBumsl 
|et al.|2004 ^, and needn't have been properties of its formation, as 
is supposed for the Galilean satellites. For a recent review of the 
issues surrounding satellite formation in circumplanetary discs, see 
[Estrada et alT|p008l l. 

It should be noted that we do not find circumplanetary discs 
around protoplanets with masses similar to Uranus and Neptune. 
This may not be a problem for Neptune, since Triton, Neptune's 
primary satellite has a retrograde orbit. However, Uranus' satel- 
lite system exhibits properties that are suggestive of a disc origin. 
Again, this might be resolved by considering the eventual cooling 
of such a protoplanet's rotating envelope following the dispersal of 



the encompassing circumstellar disc. As the envelope cools it may 
well flatten into a disc suitable for satellite growth. 



5 CONCLUSIONS 

We have investigated the effects of grain opacity and protoplane- 
tary mass upon the properties of circumplanetary discs. We per- 
formed calculations for protoplanets ranging in mass from 10-333 
M0, with grain opacities from 0.1% to full interstellar values. Our 
protoplanets have radii of 1% their respective i?Hiii> excepting one 
comparison calculation of a Jupiter radius. 

As was stated in | Ayliffe & Bate[ ( |2009) we find that 100 
protoplanets at reduced opacities (10, 1, and 0.1% interstellar grain 
opacities) possess circumplanetary discs, as do all the higher proto- 
planet masses that we model at all opacities. The radial extent of all 
these circumplanetary discs are in good agreement with the analyti- 
cal prediction for discs formed in gaps evacuated by high mass pro- 
toplanets of one third of the Hill radius (Quillen & Trilling 19981. 
Such large discs might contribute significant flux to direct imaging 
detections of young exoplanets by reflecting light from the central 
star. Indeed, it has been proposed that the recent imaging of the 
planet orbiting Formalhaut (K alas et al.|2008| may have included 
significant flux from a circumplanetary disc. 

The discs produced are all thick in comparison with their par- 
ent circumstellar discs, and their scaleheights are found to vary with 
opacity. They become thicker at higher grain opacities, as the ther- 
mal pressure supports them more against the vertical gravitational 
pull of the protoplanet. This gravitational influence in turn leads to 
a scaleheight trend with protoplanet mass, the more massive proto- 
planets possessing more flattened discs. 

We have found that the temperatures achieved in the circum- 
planetary discs are high when compared with the condensation tem- 
perature of ice. Only at the outer edges of our discs, and then only in 
reduced grain opacity cases, does the temperature fall low enough 
for ice formation. This poses a problem vis-a-vis the formation of 
icy satellites around gas giant planets. 

Finally, in contrast to the three-dimensional models of circum- 
planetary discs of [Machidaj ( [2008| , we fail to find a ring-like den- 
sity enhancement within the circumplanetary disc. Rather, we find 
that the density decreases monotonically with radius from the pro- 
toplanet. 
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